
/%3- / JlGa f 



This Research Was Sponsored by 
THE NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


Digital Servo Loop Gain Optimization 
for 

Second and Third Order Systems 


J?DC Report -1 t-62<v10~*| 



by 

Roy Waas 


Harry W. Mergler 
Associate Professor of Engineering 
Principal Investigator 
Ns G 36-60 


Digital 

Systems 

Laboratory 



ABSTRACT 


j'vLoi 

In this paper, an optimum loop- gain -frequency function 
is derived for a continuous linear servo system of second or third 
order operating with incremental data. 

A ramp system input is closely approximated by integrat- 
ing a pulse sequence and using the resulting sequential step func- 
tion. The difference between the ramp and the system output is 
taken as error. The error squared integral is used as a criterion 
of the goodness of approximation of the ramp by the system output 
and this integral is minimized for second and third order linear 
systems. Digital computer techniques are used to minimize this 
integral. 

In both the second and third order systems, the functional 
relationship between the loop gain and frequency which minimizes 
the error squared integral is found to be linear. Furthermore, 
the family of linear minimizing functions generated in the third 
order case is of the same slope as the second order minimizing 
function. It is also shown that by choosing the loop gain and fre- 
quency correctly, it is possible to obtain values of the error 
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squared integral very close to zero provided the pulse sequence 
is not delayed more than a small fraction of the pulse period. 
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INTRODUCTION 


The purpose of the following analysis is to derive an 
optimum functional relationship between the d. c. loop gain and 
pulse frequency for a linear continuous servo system using incre- 
mental data. Second and third order systems operating with a 
sequential step input will be analyzed. 

A block diagram for a general linear feedback system 
is shown in Figure 1. The transfer function for this system is: 


g(» ) _ G(s) 

R(s) ~ 1 + G(s) H(s) 


( 1 - 1 ) 


In the following analysis - 

I. System output 0(t) is required to describe, as closely as pos- 
sible, a linear function of time. The desired output is: 


Y(t) = ct. (1-2) 

II. The reference input r(t) is a sequential step function (inte- 
grated pulse train). See Figure 2. 
k 

r (t) = Sp(t-kT) (1-3) 

0 
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III. Actual output 6 (t) will also consist of a series of steps. How- 
ever, their shape will be determined by physical characteristics 
of a particular system. 

IV. The first unit step is applied with the system at rest; but, in 
general, succeeding steps are applied with variable initial con- 
ditions which depend upon the value and slope of 6 (t) at the instant 
the step is applied. See Figure 3. 

V. The difference between the actual output d (t) and the desired 
output Y(t) will be called error. 

e (t) = 0 (t) - Y(t) (1-4) 

As a criterion of the goodness of approximation of the desired 
output by the actual output, the error squared integral 



is used. K is taken large enough to assure steady-state behav- 
ior; steady- state in the sense that initial conditions will not vary 
essentially from pulee K to pulse K + 1. Reasons for this choice 
of criterion are discussed in Sokolnikoff; "Advanced Calculus", 


page 378. 
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In order to simplify the analysis, the first unit step input 
occurs at t = 0, Transfer functions have also been selected in a 
form that will simplify equations. 

VI. The problem is to find a loop-gain-frequency relationship 
which minimizes I. 



DERIVATION OF THE ERROR SQUARED INTEGRAL 
FOR THE SECOND ORDER SYSTEM 


The transfer function for a second order underdamped 
system can be written in the form: 


0(s) 



(s + a) 2 + p 2 


( 2 - 1 ) 


2 2 2 

where a and p are constants and p Q = a + P • R(s) can be 
found by transforming equation 1-3 into the domain of the Laplace 
variable s. 

R(s) = — [1 + e“ St + e" 8t (2-2) 

s 

Substituting equation 2-2 into equation 2-1: 
p 2 

e(s)= “5 2 • "T [l + e" St + e" 2st + . . . . .] (2-3) 

(s + a; + p Z 


Transforming equation 2-3 into the time domain: 

P 0 _„t 

® (t) = 1 + -p— e sin[ pt - <j> ] 

+ 1 + e a (t T) p (t-T) - <j> ] 


4 


5 


+ 


+ 1 + -jj- e a ^ KT ^ sin [ p (t-KT) - <j> ] 

Pq K . % 

0(t) = K + l+— S e “ a(t ' kT > gin [p(t-kT) - <fi] , 


P 


k=0 


(2-4) 


-1 p 


whe re <j) = tan 

-a 


Changing the sine function to exponential form, 
e -a(t-kT) g . n [p( t _ kT ) . ^ j 

= e " at . e akT [e jt P (t " kT) ~ *1 . e " j [ P (t " kT) " * j 

2j 

= e’ at .e j[pt " ^^. e akT .e ‘ jPkT 
2j 

_ e -at >e j[pt-^] >e akT e jpT 


A summation of terms involving k from equation 2-5 can be ex- 
pressed in closed form since they form a geometric series. 


s e kT(a + jp) =1 _ e (K + l) (a +jp)T 
c=0 x _ e (a + jp)T 


( 2 - 6 ) 
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2 .»(*-») 

k=0 


l e (K + l) (a -j{3)T 

i_ e <“ - JP> T 


(2-7) 


At this point let a new time variable x be assigned with zero 
reference at t = KT. Then 


t = KT + x. (2-8) 

Making use of equations 2-6, 2-7, and 2-8, the summation of 
equation 2-4 becomes: 


2 e“ a (t “ kT > sin [ p (t-kT) - <f> ] 
k=0 


1 

2j 


’ e -a(KT + x) <e j[p(KT + x) -* ] 1<? (K + l) (a - j§)T 

l-e (a “^> T 


_ e -a(KT + x) >e -j[p(KT + x) - *] l e (K + l) (a + jp)T * 

1 _ e (a+jp)T p 


2j 


e -KT(u -jp). e -x(a-jp) e -j^ l e (K + l)(a -jp)T 


_ e - KT ( a + jp). e -x( a + jp) e j^ 


1-e (K + l)--(a + jp)T ' 

i . e ( a+ ) p) t . 


(2-9) 
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Expressing <jt in terms of a and p. 


= cos [tan ] + j sin [tan ] 

-a -a 





[a-jp] . 


It follows that 


( 2 - 10 ) 


e 


~3<t> _ 



[a + jp ] . 


( 2 - 11 ) 


Using equations 2-10 and 2-11, equation 2-9 becomes 


S e" a ( t-kT ) S in [p(t-kT) - 0 ] 
k=0 


a + jP „-(a -jp)x e - KT(a-jP) _ e ^ a- jP) T 

Z &0 6 Ue (a -JP)T 


+ 


a- JP 
2 &0 


e -(a + jp)x 


e -KT(a + jp) _ e (a + jp)T 

l- e ( a + ^ T 


( 2 - 12 ) 


Substituting equation 2-12 into 
assure steady state response, 

e (x) = K + 1+ iiiL 
AJP 


equation 2-4 and taking K large to 
the output becomes 
e -(d- jp)x 

• e -(«-jf»)T 
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+ 


Q_-j_P 


_ e -(a + jp)x 

e- (o + j P > T -l 


(2-13) 

The slope of the desired output function y(t) is ~ ■ . From equa- 
tions 1-2 and 2-8, it follows that 


Y(x) = K+ 


(2-14) 


Subtracting equation 2-14 from equation 2-13 and squaring the 
result, the error squared integral for the second order system 
becomes: 



2x 

T 


+ (a + iP) 2 

-4|3 2 


e -2(a- jp)x 


t - )P > 2 

-4P 2 


e ~2(a + jp)x P 0 2 

( e -<“ + i « T -,) 2 + I ? 


-2ax 


(e -(«-JP)T_l) ( e -(a. + JP)X_ 1) 


f 1 - — 


|7(«+JP) 

u- jp) .- <<,+j(5)x 1’ 

l Tj 


_\ jP(e” (a ' j( 3 )T -l) / 

, jp(*~ a+ 2 ^^-i) j\ 


r T e- ( -«»*dx=- 
Jo ( tt “jP) 





(2-15) 


(2-3#) 
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x -(a - j(3)x 

a-jP 


T 

0 


+ 


1 

(a-jP) 



T e “( a ~ jP) x 


dx 


T 

a- jp 


e -(a-jP)t _ 1 ( e ’ (a - jP)T -l) 

(a-JP) 2 


(2-17) 


Integrating equation 2-15 with the aid of equations 2-16 and 2-17: 


T 3 2T 2 

I = 1+ — ~~ 2 . + 

3T 


(q +JP) 


2T“ jpT(e -(q^jP) T _i) 


(! - e -( a -^ T ) 

(a - JP) 


S an, ig) 


(l.e-<°*fr> T ) 


jpT (e-< a + ^) T - 1 ) ' 


+ (a-t-JP) (-D 

jPT 2 (e~ ° “ j ^ T -l) 


Te 


-(a-jP)T 


a- jp 


(«-jP) 2 


+ 



(q- jp) 

(e-<‘ + i« T -l) 


a +jp 


(e -(. + jP)T_ 1 ^ ~ 
(• + JP) 2 


+ 


8p 2 T(a -jp) (e ^ a " jP ^ T -l) 2 


+ 


(a- JP) 2 (e~ 2(a + j p)T -l) 
8p 2 T(a +jp) (e‘^ a +jP)T -i) 2 



3 iPT( “- jW jpT(o+jW ' iPxVjp) 2 JOtW 2 


jPT(a-jp) (e' ( “' J|i)T -l) jPT(a + jp) (e‘ (a + jP)T -l) 


(a + jP) 2 (e~ 2 ^ tt ~j^ T -l) 


8P T(a-jp) (e 


-(a-jp)T 772 


(a-JP ) 2 (e- 2 ;^^ T -l) 

8p 2 T(a + jp) (e ~^ a + j ^ T -l) 2 


P 0 (e" 2aT -l) 


4a p 2 T(e-(a-jP) T _i) (e -MjP)T_ 1} 


(2-19) 


Equation 2-19 can be simplified by combining terms which appear 
in conjugate pairs. 


2 R e a - 2 Re ( a + *P) 

-jPT(a-JP) " 2Re 2 

. -JPTP 


4a Pj 
■JPTP n 


( 2 - 20 ) 
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2 Re 


n + jP 

j|3T 2 (a-jf3) 2 


= 2 Re 

(a + jp) 3 
2 4 


JP T 0 O J 


= 2 


3a Z p-p 3 

0 T V 


, 2 2 
6a -2p 

R 4 T 2 

Pq T 


( 2 - 21 ) 


2 Re 


(„ + jP) 


LjPT(a-jP) (e' ( “' J|3)T -l) 


= 2 Re 


(a +M. 


LjPTP 0 2 (l-e (a_jP)T ) 


= 2 Re 


-j(a 2 -p 2 + 2gpj) 


pp 2 T(l-e aT cos pT-je aT sinpT 




2 Re 


(n.+ jP) Z [(-j) (l-e aT cos pT) - e aT sin pT] 
PP q 2 T [ (1 - e aT cos PT) 2 + e 2aT sin 2 pT] 


2 2 a_T T* 

-(a -p ) e sin pT-2gp (e a cos PT - 1) 

PP 0 2 T [(l-e aT cos PT) 2 + e 2aT sin 2 pT] 


2 2 n.T nr 

2(p -g ) e sin pT + 4aP (1- e a cos pT) 

PP q 2 T [l+e 2aT -2e aT cos pT] 


( 2 - 22 ) 
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2 Re 


(g + jg) Z (e~ 2(a ' jP)T -l) 
8P 2 T (a-jp) (e‘ (a ' j|5)T -l) 2 


= 2 Re 


(a + jp) 3 [ 1 + e a '^' cos PT + je a ^sin pT 
[ 8p 2 p Q 2 T [e" aT cos PT + je" a T sin p T - 1] 


= 2 Re 


(a + jP) 3 [e 2a ^'cos 2 pT - 1 + e~ 2a ^’sin 2 pT - 2je a ^sinpT] 


8p 2 p Q 2 T [(e“ aT cos PT-1) 2 + e" 2aT sin 2 pTj 


= 2 Re 


[a ^-3gp 2 + j(3g 2 p ~p 2 )] [e 2aT -l-2je aT sin PT] 


8p 2 p Q 2 T [e' 2aT -2e" aT cos PT + 1] 


(a 3 - 3gp 2 ) (e~ 2aT -l) + 2(3g 2 p - p 3 ) e° aT sin pT 
4p 2 p Q 2 (e' 2aT -2e" aT cos pT + 1) 


(2-23) 


P 0 (e- 2 aT -l) 


4 a p 2 T (e“ (a “ jP)T -l) ( e ~ (a + jf3 > T -l) 


PO (l-e- 2 * T ) 


4 a p 2 T [e" aT (cos pT + j sin pT) -1] [e“ a T (cos pT- j sin PT)-1] 


^0 (1 - e~ 2a T ) 


4 a p 2 T [ (e~ aT cos PT - l) 2 + e“ 2aT sin 2 pT] 
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Po Z (l-a- 2 ° T ) 

4ap 2 T [e" 2aT -2e" aT cos pT + 1] 

(2-24) 

Substituting the results of equations 2-20 through 2-24 into equa- 
tion 2-19: 



PP Q 2 T (l + e 2aT -2e aT cos pT) 


+ (a 3 - 3qp 2 ) (e~ 2aT -l) 4 2(3q 2 p - p 3 ) e~ aT sinpT 

4p 2 p Q 2 T (e‘ 2aT -2e" aT cos pT + l) 

+ P 0 2 (1 ~ e~ 2aT ) 

4 a p 2 T (e“ 2aT -2e‘° T cos pT+l) 

J_ 4q 6a 2 - 2p 2 + 2(p Z -q 2 ) sin pT+ 4qp(e~ aT -cos pT) 
3 P q 2 T P q 4 T 2 2PP q 2 T (cosh aT - cos PT) 

, (3 a p 2 - a 3 ) 2 sinh aT + 2(3 a 2 p - p 3 ) sin pT 

+ 2 2 

8p Pq • T(coshaT - cos pT) 
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2f3g sinh aT 


8a P T (cosh aT - cos pT) 


Simplifying, 


(2-25) 


I 


1 _ 

3 


. ^ 2 oa 2 

4a 6a - 2p 

_ 2^ Q 4_2 
|3 0 T T 


, (P + 5a p) sinh aT + (3a P - a ) sin pT + 8a P (e -cos pT) 

2 

4a pp T (cosh aT - cos pT) 


( 2-261 


Equation 2-26 expresses the error squared integral I as a 
function of three variables; a, p, and T. It can now conven- 
iently be transformed into a function of the scaled time variable 
aT and damping ratio t, . The two forms of system character- 
istic equation can be compared below. 

S 2 + 2£W Q S + W Q Z = 0 
S 2 + 2aS + a 2 + p 2 = 0 
It can be seen that 


2 2 2 2 2 2 
W Q = a + P and a = i W . 
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Then 


2 2 , 2* 2 1 - t 

P = W Q (1 - & ) = a — 2^- 


Let 


i-r 


= A 


(2-27) 


If the numerator and denominator of each term of equation 2-26 
are multiplied by the appropriate power of T and the substitutions 


aT = Z and £T = AZ 


(2-28) 


are made, I takes the form: 


1 


4Z£ 2 (6Z 2 - 2A 2 Z 2 ) t, 4 

1 " 3 " 2 + 4 

Z Z 


(A 3 Z 3 + 5AZ 3 ) sinh Z 4 (3A 2 Z 3 - Z 3 ) sin AZ + 8AZ 3 (e~ Z -cos AZ) 
4AZ 4 (1 + A 2 ) (cosh Z - cos AZ) 


_L ii! + 2^ 4 (3 - A 2 ) 
3 ' Z ^ 
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+ (A^ + 5A) sinh Z + (3A^ - 1) sin AZ + 8A (e~ Z -cos AZ ) 

4AZ , . 

— 5 — (cosh Z - cos AZ) 


From equation 2-27 


1 + A 2 = A: 


therefore 


3 - A 2 = 4 - (1 + A 2 ) = & ~ 1 , 

r 


A 3 + 5A = A(l + A 2 + 4) = A * 1 , 

r 


(2-29) 


and 


3A 2 - 1 = 3(A 2 + 1) -4 = p . 

r 

Finally, using these relationships in equation 2-29, 

i . j_ . sf . idioH - » 

' 3 z z 2 

+ (4; 2 -H) sinh Z + (3-4£ 2 ) (sinAZ) /At 8£ 2 (e' Z -co» AZ) 

4Z (cosh Z - cos AZ) 


(2-30) 


DERIVATION OF THE ERROR SQUARED INTEGRAL 
FOR THE THIRD ORDER SYSTEM 

The transfer function for a third order underdamped 
system can be written in the form: 


^(S) 7f3 0 

R(S) (S + T )t(S+ a) 2 + p 2 ] 


(3-1) 


2 2 2 

where a, |3 , and 7 are constants, | 3 Q = a + |3 , and R(S) is 

defined by equation 2-2. 


Then 


0j(S) = 


TP, 


(S + 7) [(S + a) 2 + p 2 ] S 


1 n x -ST. ~2ST 
-5- [1 + e + e + 


and 


\ 1 


+ 1- % e -T(t-T) + 2fo e -a(t-T) s . 

7 2 % 


sin[p(t - T) -<^ 1 ] 
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L . P 0 -y (t - KT) 7P 0 -a(t-KT) 

+ 1 - — - e ' + - — e 'si 

2 


sin[p(t-KT) ■‘<j> L ] 


(3-2) 


S.W-K + l -X j=„ e-<‘- kT > 

1 2 k=0 By „ k=0 


sin[ p (t - kT) , 


, ,2 „2 

(a - y) + p and <p. 


-IS -IB 

tan ■■■ — - + tan — — 

-a y - a 


(3-3) 


S .-* t - kT > s e™ 


e"^ (i . e (K + l)rT\ 


(3-4) 


From equations 2-5 through 2-9 it can be seen that 


S e a ^ sin[ p (t - kT) - $ . ] 

k=0 
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1 

2j 


-j* -KT(a-JP) -x(a-iP) 1 - e (K + *> (a " JP > T 


- e 


j^l , e -KT(o+jP) e -x(a+jp) l e (K + l) (o+jfJ)T 

l-e<“ + ^> T 


(3-5) 


where x is defined by equation 2-8. 


Let tan * — = a, and tan — 

-a 1 


7- a 


= Pj. Then from equation 


3-3 


e^l = cos (ftj + Pj) + j sin( ai + P^ 


= cos cos p^ - sin sin p^ + j(sin cos p^ + cos sinp^) 

= . _s_ ii -») _ JL _P + i 


0 Y 0 P 


0 7 0 


f P (7 - a) . j 

f-a pV 

0 

■ 

0 

CQ_ 

1 

l"o ' ^o)J 


Pq y o 


■a (7 -a) - p +j(P[7-a] - aP ) 


(3-6) 


and 


e . 1 = 


P 0 7 0 


-a (7 - a) - P - j(p[ 7 - a] 




(3-7) 
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In terms of x, the right hand side of equation 3-4 becomes 


e -Y(KT + x) (l _ e (K + l)YT ) e -Y(KT + x)_ e 7 (T-x) 


1 - e 


7 T 


1 - e 


7 T 


(3-8) 


By using equations 3-4 through 3-8 and taking K large in equations 
3-5 and 3-8, equation 3-3 can now be transformed into a steady 
state response function for the third order system. 


p 


0 


e l (x) = K+i - 


-e 


7 (T - x)' 


1 - e 


YT 


2jP7 f 


P 0 T 0| 


-a(T-a)-p 2 -j(p[T-a] -«(() e’**' 1 -'' 15 ' 


_ e (a -j(3)T 


-a (7 


a) “ P 2 + j(P[ T-a] -ap) 

j 


e “ x (a + jp) 


_ e (a+jP)T 


(3-9) 


Let 

a (7 - a) + p 2 = M , 2 a p - pY = N, 
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M - jN 

e “(a -jP)T 


= D 


and 


M + jN 
e -(a+jp)T 


D 


(3-10) 


Using equations 3-10 and 2-14, the error squared integral for 
the third order system can be expressed as follows: 



(3-11) 


Squaring the integrand: 



4 *V 


-,o 2 -7X 

2 * 2 ^0 8 
T 2. - 7 T 


Zy 


V (e 


i) 


2jp7 f 


De -*(a-jP) _ 5 e -x(a+jp) 
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_ _ 2 -7x 

2 P 0 * e __ 

T7„ 2 (e' 7T -l) jTp7 0 2 


yx 


De -x(a -jP) _ 5e -x(a+jp) 


_ 2 -yx 

Pq ^ 

jP 7 0 4 (e- 1 ' T - 1 ) 


De -x(a -jP) _ 5e -x(o + jp) 


dx 


(3-12) 


Integrating equation 3-12 with the aid of equations 2-16 and 2-17: 


*i ■ 1 + T - 


l S>o 4 < e ' 2TT - 1 ) 


4 . -7 T ,.2 , mo 2 4 
2T7 Y q (e -1) 4T£ 7 Q 


f 


d 2 - 2 (a - jP)T dB _ 2aT B 2 - 2 (. + jp)T 

— wm r (e - 1)+ — Jfrrp — 


-2p 

- 1 + ^ + — 1 


TYTq 2 jTp7 Q 2 


D( e- (a -jP )T -l) 5 (e~ (a +j P >T -l) ' 

- (a * jP) (<* + jP) 


2 P, 


1 2 7 0 2 (e' TT -l) 


i e^ T + -L (e-^-D 


.2 2 

J T P7„ 


D 


f T (e-( a -jP> T ) + (e -(«- iP ) T - i ) 


a - jp 


(a “ jP) 
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- D 


Te 


-(a+j{J)T e -(a+jp)T_ 1 


a + jp 


(a + JP ) 2 


/ J 


0 2 


jTPT 0 4 (e“ TT -l) 


^-(g+y-jP)!^ 

-(a + 7 - jp) a+y+jp 


(3-13) 


Equation 3-13 can be simplified by combining terms which appear 
in conjugate pairs. 


p , 


P Q 4 (e" TT +l) (e"' VT -l) 

j 3yT 2 ” 

2Tyy n (e 7 -1) 2Tyy, 


4 cothJ T 


(3-14) 


4Tp 2 y 0 4 


. 2 Re 


(M 2 - N 2 - 2jMN) (e“ 2(a ~ jP)T -l) 

L 2(. -») (e- (a ' iP)T -l) 2 


Arno Z 4 
4Tp y Q 


Re 


(M Z - N 2 - 2jMN) (e~ aT cos pT + 1 + je~ aT sinpT) 
(a - jp) (e" aT cos pT - 1 + je " aT sinpT) 


_ 7 2 Re f[a(M 2 -N 2 ) + 2p MN + j(p(M 2 - N 2 ) - ZaMN) 1 
4TP 2 7 q 4 P 0 2 ((e" aT cospT - l) 2 +e" 2aT sin 2 pT ] 

Continuation of numerator 

* [e~ 2aT cos 2 pT - 1 - 2je ~ aT sinpT + e‘ 2aT sin 2 pT 
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7 2 Re ([a (M 2 - N 2 ) -I- Z(3 MN + j(p (M 2 - N 2 ) -2g MN) ] 
4Tp 2 P Q 2 7 0 4 [e' 2aT -2e " aT cos pT + l] 


* [(e~ 2aT -l) -2je ~ aT sin pT j } 


f[a(M 2 - N 2 ) + 2pMN] [-sinhaT ] + sinpT [p(M 2 - N 2 ) - 2aMN]^ 

O 7 A ~ ' 


2 2 4 

4Tp Pq 7q (cosh oT - cos pT) 


(3-15) 


2 — . -2aT .. 2 #-2 _ T 2, , -2aT . 

7 DD (e -1) _ 7 (M 4 N ) (e -1) 

4Tap 2 7 Q 4 4Tap 2 7 0 4 (e‘ aT cos PT - l) 2 +e" 2aT sin 2 pT 

7 2 (M 2 + N 2 ) sinh aT 

4Tap 2 7 Q 4 (cosh aT - cos pT) 

(3-16) 


TP7, 


2 Im < 


D(e~( a - »> T .l) 


-(a-jp) ) 


7 

TP7 q 2 


f (M-jN) (e' (a ~ jWT -1) \ 
I -( a -jP)(e- ( °- jP)T -l)J 


. 2 Im 
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2g (3 M - N(g 2 - p 2 ) 


P 


0 


2t 


2 2 2 
T P 7 0 P 0 


T[ (p M - gN) (1 - e a ^cos (3T - (g M+ pr^)e ar ^sinpT] 
2e a ^' (cosh gT - cos PT) 


+ 2gftM-N(g 2 - p 2 ) 

(3-19) 

D(e ~( a * 7 ~ jP)T_ 1} " 

-(a +7 - jp) 

(M - jN)(e~( a * 7 > T cos pT - 1 + je~* a * ^siiftT) 
-(g + 7 - j(3) (e aT cos |3T - 1 + je aT sinpT) 

* 

[M(g + 7) + Nft + j[pM-N(g + 7)]] [ e ~( 2a + 7 ) T 
[(g + 7) 2 + p 2 ] [(e a ^'cos pT ~ l) 2 + e 2ar ^sin 2 pT 


-^o 7 


TP7 0 4 (e~ TT -l) 


Im 


2P 0 * 


Tp7 n 4 (e“ 7T -l) 


Im 


Pq l 

TP7 4 (e" TT -l) 


2 Im 


* +1 -(l + e~ 7T ) e" aT cos pT + j(l-e~ 7T ) e~ aT sinPTl 
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TPT Q 4 (e 



[M(g^y) + Np] [(l-e~ 7T ) e~ aT sin pT] + 
[(a + Y) 2 +p 2 ] 2e" aT [ cosh aT - cos pT] 


[pM-N (a , + y) ] [e~^ 2a be~ aT (l +e~ TT ) cos pT] 



* 

[M(a + y) + Np] [-sin pT] + [pM-N(g+Y)] 
[(g+Y)^+P^] [cosh aT - cos PT] 


yt 


[cos PT coth 2 + (e 


aT \ // -TT i\ 

£ ) /(e -1) 


(3-20) 


Substituting the results of equations 3-14 through 3-20 into 
equation 3-13: 


I l = 


3 


P, 


4 


2TYY r 


YT 

coth — + 


._.2 2 4 

m p 0 7 0 


[P(M 2 -N 2 ) -2aMN] gin gT -[q(M 2 -N 2 ) +2PMN] sinh qT 

cosh gT - cos pT 


+ 


4Tgp 


Y 2 (M 2 + N 2 ) sinh gT 
2 4 , 

Y q (cosh gT - cospT) 


2P 0 2y (g N - pM) 

2 2 2 
T ^0 T ^0 P 0 
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ifo_ 

T V 


T + 1 ‘ 

+ - 27 

7(1 - e TT ) 7 2 

2 2 

T w 0 7 o L 


2gpM-N(g 2 - p 2 ) 


■aT 


T[(PM - gN) (e -cos PT) - (aM + pN) sin PT 
2(cosh gT - cos pT) 


7P, 


Tp7, 


- [M(g + 7) + PN] sin pT + [pM - N (g +7)] 
[(g + T ) 2 + p 2 ] [cosh gT - cos pT] 


[cos pT coth \ + (e~ aT + e^ a ' l ' T ^ T ) /(l- e YT )] 


(3-21) 


Let gT = Z and PT = AZ as in the case of the second order system 
where A is defined by equation 2-27. In addition let 


yT = XZ 

Then from equations 3-1, 3-3, and 3-10: 

MT 2 = gT(y T- gT) + p 2 T 2 = Z 2 (X -1 + A 2 ), 

NT 2 = 2gpT 2 - PtT 2 = Z 2 (2- X) 

T q 2 T 2 = (gT - tT) 2 + p 2 T 2 = Z 2 [(1 - X) 2 + A 2 ] , 


(3-22) 


29 


and 


o 2 2 

P„ T 


2 t 2 4. 
o T + 


2 2 
P T = 


z 2 (i 


A 2 > = 


z 2 /; 2 


(3-23) 


Using equations 2-27, 2-28, 3-22, and 3-23, equation 3-21 can 
be expressed as a function of £ , X , and Z. 


P 0 A 7T Z 4 /r coth — 

— coth — = ^ - 

2yy T 2XZ-Z [(l-\) +A ] 


coth 


XZ 

2 


2 \Z£ 4 [ (1- X) 2 + A 2 ] 2 


(3-24) 


7 2 [ P(M 2 - N 2 ) - 2q MN] 

^\ 2 y 0 4 


* 

X 2 Z 2 [AZ(Z 4 [X-1 + A 2 ] 2 -A 2 Z 4 [2- X] 2 ) -2Z 3 [X-1 + A 2 ] 

22 2 4 r . ,2 2,2, ,2. 

4A Z * Z • Z [(1- X) + A ] (l/£ ) 


■az 2 [2-x]1 
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?;Vr(\-l + A Z ) 2 - A 2 ( 2 - x) 2 - 2 (X -1 + A 2 ) (2-X) j 

4AZ [ (1 - X) 2 + A 2 ] 2 


(3-25) 


-7 2 [g (M 2 - N 2 ) + 2pMN 
._ 0 2. 2 4 

4T P P 0 Y 0 

\ 2 Z 2 [Z(Z 4 [\ -1 + A 2 ] 2 -A 2 Z 4 [2- x] 2 ) + 2AZ 3 (X-1+ A 2 ) 
-4A 2 Z 2 (Z 2 /; 2 ) Z 4 [(l-\) 2 + A 2 ] 2 

* -AZ 2 [2 - A] ] 


^ 2 X 2 [(\-1+ A 2 ) 2 -A 2 (2 - X) 2 + 2A 2 (2 -X) (X- 1 + A 2 )] 
4A 2 Z [(1-X) 2 + A 2 ] 2 


(3-26) 


7 2 (M 2 + N 2 ) _ X 2 Z 2 [Z 4 (X- 1 + A 2 ) 2 +A 2 Z 4 (2- X) 2 ] 

4aTp 2 7 Q 4 4Z-A 2 Z 2 Z 4 [(l - X) 2 + A 2 ] 2 


X 2 f(X- 1 + A 2 f 4- A 2 (2 - X) 2 ] 
4A 2 Z[(1- X) 2 + A 2 ] 2 


(3-27) 
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2 P, 


T?7, 


2Z 2 / l 2 


-2 


\Z 3 [(1 - X ) 2 +A 2 ] 


XC 2 Z [(1 - X ) 2 + A 2 ] 


2t (gN - PM) 
2 2 
p 0 


2 XZ [AZ 3 (2 -X) -AZ 3 (X- 1 + A 2 )] 

az 3 (z 2 /; 2 ) [(i-x ) 2 + a 2 ] 


2 X£ 2 [2-2X + 1-A 2 ] = 2 X £ 2 [3 - 2X - A 2 ] 

z[(i-x ) 2 + a 2 ] z[(i-x ) 2 + a 2 ] 


2 Pr 


2 2 
T 7 0 


I_ + J_ 

7(1 - e 7 T ) 7 2 


2Z 2 /l Z 


Z 2 [(1 - X ) 2 + A 2 ] 


1 


XZ(l-e XZ ) X 2 Z 2 


27 [2gpM - N (a 2 - P 2 ) 1 

T 2 pp 2 7 2 p 2 
HP 0 T 0 p 0 


(3-28) 


(3-29) 


2 ^ , ,2 2 , 

C XZ [(1 - X) + A ] 


1 - e 


XZ 


1 

XZ 


(3-30) 
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2\Z [2Z AZ -Z 2 [X - 1+ A 2 ] - AZ 2 (2 - X) (Z 2 - A 2 Z 2 )1 
AZ (Z 4 / C 4 ) Z 2 [(1 - \) 2 + A 2 ] 


2 & 4 X [2(\ - 1 + A 2 ) - (2 - X) (1 - A 2 )] 
Z 2 [(1 - X) 2 + A 2 ] 


(3-31) 


2?T (pM - aN ) = XZ [AZ 3 (X- 1 + A 2 ) - AZ 3 (2 - X)] 

2T 2 PP 0 2 t 0 2 az (z 2 /; 2 ) Z 2 [(1 - x ) 2 + A 2 ] 

XC 2 A [2X- 3 + a 2 ] 

AZ[(1 - X) 2 + A 2 ] 

(3-32) 

-2yT (gM + gN) = -XZ [Z 3 (X- 1 + A 2 ) + A 2 Z 3 (2- X)] 

2x2 PP 0 2 T q 2 AZ(Z 2 / ; 2 ) Z 2 [(1 - X) 2 + A 2 ] 

£ Z XrX- * 1 + 3A 2 - 2A 2 X] 

AZ [il -X) 2 + A 2 ] 

(3-33) 

T ^0 M(a + t) + NP 

tPt 0 4 (o + t) 2 + P 2 
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xz(z 2 /c 2 ) [Z 2 (X - 1 + A 2 ) (Z + XZ) + a 2 z 3 (2 - X)j 
AZ 5 [(1 - X) 2 + A 2 ] 2 [(Z + XZ) 2 +A 2 Z 2 ] 


XZ [ (X - 1 + A 2 ) (1+ X) + A 2 (2 - X) ] 
; 2 A [(1 - X) 2 + A 2 ] 2 Z 2 [(l+ X) 2 +A 2 ] 


(3-34) 


-7P 0 2 [PM - N (a + t)] 

TP7 q 4 [la + 7) 2 + P 2 ] 

-XZ (Z 2 / £ 2 ) [AZ 3 (X - 14 A 2 ) - AZ 2 (2 - X) (Z - XZ)] 

AZ -Z 4 [(1 - X) 2 + A 2 ] 2 [(Z + XZ) 2 + A 2 Z 2 ] 

AX Z [X 2 + A 2 -3] 

A £ 2 [(1 - X) 2 + A 2 ] 2 Z 2 [(1 + X) 2 + A 2 ] 

(3-35) 

Substituting the results of equations 3-24 through 3-35 into 
equation 3-21: 


1_ coth (XZ/2) 

3 2XZ £ 4 [(1 - X) 2 + A 2 ] 2 
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* 

+ A£ Z \ 2 [(X - 1 4 A 2 ) 2 -A 2 (2 - X) 2 -2(X - 1 -f A 2 ) (2 - X)] sin \Z 
4A 2 Z [(1 - X.) 2 + A 2 ] 2 [cosh Z - cos AZ] 

* - &V [(X - 14 A 2 ) 2 -A 2 (2 - X) 2 4 2A 2 (2 - X) (X - 1 + A Z )]siA Z 

+ — t (X ~ 1 + A 2 ) 2 + A 2 (2 - X) 2 ] sinh Z 
4A Z Z [(1- X) 2 + A 2 ] 2 [cosh Z - cos AZ] 


2 

X c 2 Z [(1 - X) 2 + A 2 ] 


+ 


2X ; 2 (3-2X - A 2 ) 
Z [(1 - X) 2 + A 2 ] 


XZ£ 2 [(1 - X) 2 + A 2 ] l-e XZ XZ 


+ 2X k 4 t 2 ( X ~ 1 + A 2 ) - (2 - X) (1 - A 2 )] 

2 r , .2 2 , 

Z [(X - 1) + A ] 


+ ; Z AX(A 2 -3 + 2X) (e~ Z -cos AZ) -£ Z X (X - 1 + 3A 2 - A Z X) sin AZ 

AZ [(1 - X) 2 + A 2 ] [ cosh Z - cos AZ] 

+ X(3A 2 +X -1) sin AZ - AX [A Z + X 2 -3] [cos AZ coth (XZ/2) 

A £ 2 Z [(1 - X) 2 + A 2 ] 2 [(1 + X) 2 + A 2 ] [cosh Z - cos AZ] 


* 


+ (e 


- Z Z(1 + X) 


+ e 


USLl 


xz., 

* R 


(3-36) 
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In order to simplify equation 3-36 and to facilitate programming 
for digital computation, let 


X - 1 + A 2 * = H, 

2 - X = P, 

(1 - X) 2 + A 2 = G, and 


cosh Z - cos AZ = R. 


(3-37) 


Then 


- _ _1_ coth (XZ/2) 

1 3 2xz; 4 g 2 


+ a £ 2 X. 2 [ H 2 - A 2 P 2 - 2 PH ] sin AZ - t, 2 X 2 [ H 2 - A 2 ^ + 2A 2 PH] sinh Z 


2 2 
4AZGR 


2 r 2 22, 

X [H + A P ] sinh Z 

2 2 

4A ZG R 


X£ 2 ZG 


2X£ (P-H ) _2 

ZG xz C 2 g 


1 - e 


i + J- 

xz xz 
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. 2X£ 4 [2H-P(1- A 2 )] 

z 2 g 

AX £ 2 (H - P) (e" Z -cos AZ) - (H + A 2 P) sin AZ 

AZGR 

2 2 * 

X [3A + X - 1] sin AZ - [G - 2P] AX [cos AZ coth (XZ/2 ) 

AC 2 ZG 2 (G + 4X)R 

* Me" 2 * e Z ^" X W-e XZ )l 

(3-38) 


INTRODUCTION TO THE PROGRAMMING PROGRAM 


Equations 2-30 and 3-38 represent the error squared 
integrals for second and third order underdamped systems 
respectively. If the systems are critically damped. 


4 = 1 , 


A = 0, and lim 
A*0 


sin AZ 
A 


= lim Z cos AZ = Z 

Ar*0 

by L.' Hospital's Rule. 


(4-1) 


In this case, equation 2-30 becomes 

1 _4_ 6 5 sinh Z - Z + 8(e ^-1) 

3 “ Z z 2 4Z (cosh Z - 1) 

(4-2) 


and equation 3-38 becomes 


T _1_ x (XZ/2) , \ 2 H (H - 2P) Z + 2\ 2 p (P - H) sinh Z 

1 3 4 4 

2X.ZH 4ZH R 
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Z_ + 2\(P - H) + 


XZH 


ZH 


XZH? 


1 - e 


i— + — 

XZ \Z 


2X [2H - P] X(H - P) (e' Z -l) - \HZ 

2 2 2 
Z H ZH R 


K(X 2 - 1) Z - (\ 2 - 3)X [coth (XZ/2) + (e~ Z + e Z ^ * ^)/(l - e XZ )) 


ZH 4 R (H 2 + 4X) 


Then 


(4-3) 


If the systems are overdamped, t, >1 and A is imaginary. 


A = j 


. V, 2 .! 


= iA 0 ' 


cos jA Q Z = cosh A Q Z, and 


-j sin jA Q Z = sinh A Q Z. 


For the second order overdamped case, 

I = JL _ id + 2C 2 (4£ 2 - 1) 


(4-4) 


(1 + 4C, 2 ) sinh Z + (3-4£ 2 ) (sinh A Q Z)/A 0 + 8£ 2 (e" Z -cosh A Q Z) 

4Z (cosh Z - cosh A^Z) 


(4-5) 
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and for the third order overdamped case 


I 


1 



coth (XZ/2) 
2XZ£ 4 G 2 


+ AX 2 ^ 2 [H 2 + A 2 P 2 - 2 PH] sinh AZ + ^ 2 \ 2 [H 2 + A 2 P 2 - 2A 2 PH)sinKZ 

4ZRA 2 G 2 


, X 2 [A 2 P 2 - H 2 ] sixth Z 2 2X£ 2 (P - H) 

+ 2 2 2 + ZG 

4A ZRG XC ZG 


XZG£ 


XZ 


1 - e 


XZ 


2X^ 4 [2H - P (1 + A 2 )] 
2 

Z G 


AXC 2 [H - P] [e~ Z -cosh AZ] - ^X [H - A 2 P,jsinh AZ 

AZGR 


X j\ 2 - 3A 2 - 1] sinh AZ - [X 2 - A 2 - 3 ]AX [coshAZ coth(XZ/2 ) 
AC 2 ZRG 2 (G + 4E) 


+ (« 


+ e 


Z(1 4 X)) , (1 . e XZ„ 


(4-6) 


Consider first the second order system. The problem 
is to find a d. c. loop gain-frequency function which minimizes I, 
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or to find an equivalent relationship of the form £, = f(Z) 

which minimizes I. This can be done by setting the derivative 

of I with respect Z equal to zero and solving the resulting ex- 

2 . . 

pression explicitly for Z, = f(Z). Even the simplest of the 
error squared integrals is too complex to handle analytically 
in this manner. 

Because of the complex structure of the I functions, 
digital computer techniques were used to plot and locate the 
minimums of these functions. All programs were written in 
the language of Burrough's Algebraic Compiler. A discussion 
of the programming techniques used is taken up on the following 
pages. A listing of all programs used appears in the appendix. 



PROGRAMMING FOR THE SECOND ORDER ERROR INTEGRAL 


The evaluation of I for several values of the parameter 
t, was simply a matter of substituting values into the derived 
functions I for £>1, t, =1, £<1. Rough plots of the I family of 
curves were made first using slide rule calculations. These 
plots were used to select points for the input arrays. There 
were unexpected irregularities in some of the curves, especially 
those of very low Z, which made it necessary to rerun the three 
programs for additional points. Resulting curves are shown in 
Figures 5 and 6. 

The second phase of the problem is that of accurately 
locating the true absolute minimum for each curve. An analytic 
expression for the derivative of I with respect to Z is required 
in the minimization procedure used. The three forms of I' (Z) 
are listed below for reference. 


Underdamped I' (Z) = - 4 ^ . ~ *) 

z z z 3 


+ { [ (1 + 4Z, 2 ) cosh Z + (3-4£ 2 ) cos AZ + 8£, 2 (A sin AZ-e Z ) ] 
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Z [ cosh Z - cos AZ] - [(l + 4£ 2 ) sinh Z 
+ (3 - 4£ 2 ) ■ Sin A AZ + 8 C 2 (e" Z -cos AZ) ] [cosh Z 

+ Z sinh Z - cos AZ + AZ sin AZ]| /4Z 2 [cosh Z - cos AZ] 2 

(4-7) 

Overdamped I» (Z) = ^ ^ 

z 2 z 3 

+ {[(1 + 4£ 2 ) cosh Z + (3 - 4£ 2 ) cosh A Q Z - 8£ 2 (A q sin A Q Z 

-Z 2 

-e ) ] Z [cosh Z - cosh A Q Z] - [(1 + 4£ ) sinh Z 

sinh A Z 

+ (3 - 4C, ) + 8£ (e -cosh A.Z)] [cosh Z 

0 0 

+ Z sinh Z - cosh A q Z - A Q Z sinhA Q Z]J /4Z 2 [cosh Z - coshA Q zf 

(4-8) 

Critically damped I' (Z) = — 

z z 

+ |[5 cosh Z - 1 - 8e Z ] Z [cosh Z - 1] - [5 sinh Z 
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- Z + 8(e" Z -l)] [ cosh Z + Z sinh Z - 1] J /4Z Z [cosh Z - l] 2 

(4=9) 

A plan for accelerating the convergence of iterative processes 
described by J . H. Wegstein in Communications of the Associa - 
tion for Computing Machinery, Volume 1, number 6, page 9, of 
June, 1958, was used to find the roots of the family of curves 
I'(Z) = 0. The procedure is explained below. 

1. Write I’(Z) = 0 in the form 

Z = Z + gl* (Z) = F(Z), where g ^ 0 is a constant. 

2. Choose a Z^ (refer to Figure 7, page 59) 

3. Let Zj = F(Z q ) 

4. Let Z = F(Z ) 

^ 1 

Instead of using Z to continue the process, find the intersection 

of the secant joining Z^, F(Z^) and Z^, F(Z^) with the line y = Z 

and call this value Z . For each new Z, follow this same pro- 

cedure. Each Z _ is obtained from the intersection of y = Z 
n + 1 

and the secant joining Z , F(Z ) and Z F(Z _). Note 

& n n 7 n-1 n - 1 

from the graph that each F(Z ) = Z , . From the geometry 

n' n + 1 

of the graph, there are two relationships that are of use. 
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Z 


n + 1 


Z n + 1 Z n - 1 
Z n + 1 Z n + 1 


-Z Z 
n n 

-Z -Z 
n n 


and 


(4-10) 


Z 


n + 1 


F(Z ) 
n 


(4-11) 


A flow chart for the procedure appears in Figure 8. 

Several tries were required before the process would 
work satisfactorily. The problem was to find the right combina- 
tion of g (the arbitrary constant multiplier of I* (Z) and h (for the 
convergence test J Z^ - Z^ J <h). Too large a g caused the 
hyperbolic functions to go out of range and too small a g caused 
overflow in the calculation of Z . h also had rather narrow 
limits. Results were not reliable unless h was kept near some 
very small value, yet when h was made too small, Z oscillated 
without converging. It is also evident from the plots of l(Z) 
that it was necessary to have a good idea of where the absolute 
minimums were located since there are several roots to the 
derivative equations in some cases. Results of the minimi- 
zations are shown in Figure 9. The surprising fact is that 



45 


£ Z = f(Z) is a linear function, therefore no curve fitting problem 
exists! The equation of the minimizing function is: 


l Z = 0.256Z. 


(4-12) 


PROGRAMMING FOR THE THIRD ORDER ERROR INTEGRAL 


The third order case is complicated by an additional 
parameter X. Evaluation of 1^ was performed by using the 
same procedure that was used in the second order case except 
that data was read into the program by using the FOR statement 
rathe r: than by an input array. Several sets of points were 
needed; one set for each family of curves, since the minimum s 
form a surface rather than a plane curve. Plots of 1^ (Z) 
appear in Figures 10, 11, and 12. 

An examination of the third order error squared integral 
will reveal why it was impractical to use Wegstein's method, 
which requires the derivative of 1^ (Z), for locating the minimums 
of each function. They were found instead by constructing a net 
from the plots of 1^ (Z) in the neighborhood of each minimum. 
Spacing of the net was decreased until 3 points were determined: 
Zq, Zq + A, and Z^ - A such that 1^ (Z^) was less than either 
1^ (Zq + A) or 1^ (Zq *• A) where A is the tolerable error in Z. 
Results of this analysis are shown in the family of straight lines 
in Figure 13. Each line represents a minimum function = f(Z) 
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for a particular value of X. . These lines are parallel, therefore, 
the minimizing surface can be expressed as 

£ 2 = mZ + f(b) , 

2 

where b is the ordinant of the t, axis intercept for each line, and 

m is the slope of the family of straight lines. 

By determining a functional relationship b = f(\), an 

2 

analytical expression for t, = f(Z, X) can be written. A plot of 
b = f(X) is shown in Figure 14. The resulting curve forms one 
leg of an equilateral hyperbola. Its equation is: 

Xb = . 515 , 

(4-13) 

therefore 




2 


. 256Z + 


.515 

X 


(4-14) 


is the required minimizing function for the third order system. 

There is a very slight scattering of points in Figure 13. 
This can be attributed to the complex function which was being 
evaluated. When such a large number of computations involving 
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extremely large variations in magnitude are required to plot a 
single point, small errors are bound to occur. 


CONCLUSIONS 


t/T 

(5-1) 

t/T + b 

(5-2) 


Mathematically, a delay in the system input pulse 
sequence is equivalent to shifting the desired output y = t/T 
as shown in Figure 4. The analysis of this paper was based on 
a linear desired output (equation 5-1). Optimizing functions 
were derived based upon minimization of 



(5-3) 


It can be seen from plots of this integral that in all cases it 
approaches zero for some T. 

If a similar study is made based on equation 5-2 for a 
delayed pulse sequence, the integral to be minimized becomes: 


50 



a 


(K + 1)T 
[e(t) - t/T]‘ 


dt 




(K + 1)T b 2 

[e(t) - t/T] dt + 

T 



(K + 1)T 
dt 


(5-4) 


The first member of the right hand side of equation 5-4 has been 

plotted for second and third order systems in sections 2 and 3 of 

2 

this paper. The value of the third integral is b^ . The second 
integral can be handled with the aid of equations 2-8, 2-13, and 
3-9. In the steady-state second order case: 

2b f (K + 1)T 2b r TV 

- J [ e(t ) - t/T] dt = - — i ( il - + 

KT Jo ^ A 



e -(a-j|3)T_ 1 


(a 


.iPle^^ 


-(a + jfi)T 


-1 



dx 



T 2 a+jP 

2T 2jP[-(a-jP)] 


a - jP 

2JP f-(a + jP)] 


2b 

1 

1 2)TPP 0 2 


[(a - jp) 2 - (a + jp) 2 ] 



+ 



a 


T 


(5-5) 
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If the terms added to the integral I by consideration of a pulse 
delay are denoted by I , then for the second order case 


X „ I+ ^ 

‘ 1 1 0 O T 


(5-6) 


In terms of £ and Z, equation 5-6 becomes: 


I = b - b, + 4b, 
s 1 1 1 


Z^l + A ) 


- b 2 + 4 b n 2 

2, b l Z 


The third order case can be handled similarly. 


(5-7) 


2b f (K + 1)T 2b f Tl p 


-7X 


7 n 2 (e" TT -l) 


2jpT Q 2 


(M - jN)e~^ a ~ (M + jN)e~^ a + 


s -(a-j{3)T_ 1 


; -(a + jp)T_ 1 


dx 


2b l 

T 2 

T - 4 - 

p 0 2 

1 \ A T 

M- jN 

M+jN 


T 

2T + 

2 ' 

^0 

T 1 2 

i / 2 jPT 0 

"(a - jp) 

r ^TTJp) 
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2b S 

= - b . + -O. 

T7 0 2t 


2b 1 y(2) 

2jTpy 0 2 


Im 


M - jN 
a - jP 


M + jN 
a + jP 


(5-8) 


P, 


Im [(M - jN) (a + jp) - (M + jU)(a - jp)] 


[pM - aN] 


(5-9) 


Substituting equation 5-9 into 5-8, I for the third order case 
can be written 


2 2b A 4 V 

'l.'b, - b l + Hi - + 2—2 


TT 0 T PTP 0 Y q 


In terms of £ , Z, and X, equation 5- 10 becomes: 


(5-10) 


2 bZ 2 (l + A 2 ) 

I. =b - b. + = 5 5- 

X Z ■ Z [(1 - X) + A ] 


4b^X Z [AZ -Z 2 (X - 1 + A 2 ) - Z -AZ 2 (2 - X)] 
AZ- Z 2 (1 + A 2 ) Z 2 [(1 - X) 2 + A 2 ] 
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2b 


= b. 


V 


1 


\t, 2 z [l/t, 2 + X 2 - 2\] 


^XC 2 [^2 + 2X -4] 
Z [-^ + X 2 - 2X] 

r 


(5-U) 


Plots for equation 5-7 are shown in Figures 15 through 20 for 
several values of £ and for b ^ = 0.25, 0.50, 0.75. A summary 
of the results of varying are shown in Figures 21 and 22. 
Increasing beyond 0.25 makes it impossible to obtain a zero 
value for I. This effect is even more pronounced as is in- 
creased. Even for Z, =0.5, no clearly defined minimum exists 
for b^>0. 5. In all cases where it is possible to obtain high 
accuracy (I = 0), the optimum frequencies are lowered consider- 
ably by the pulse delay. Pulse input delays do, however, flatten 
the I curves in the vicinity of the minimums so that when I can 
be made to approach zero, gain-frequency selection is less 
critical when delay is involved. 

Equation 5-11 is plotted in Figure 23 for \ = 10, t, = 0.1 
and b^ = 0. 25, 0. 5, 0. 75. These curves should be added to the 
l , =0.1 curve of Figure 11. A casual inspection will show that 
considerable error is involved for even the b^ = 0.25 curve. In- 
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creasing either £ or makes the situation worse so that in the 
third order system not even = .25 can be tolerated. 

Plots of the error squared integral functions show that 
a linear system output can be generated to any desired degree of 
accuracy by selecting system parameters according to equations 
4-12 and 4-14. They also indicate relative errors which will 
result from not selecting parameters correctly. The slope of all 
curves increases very rapidly to the left of the minimums 
(higher frequencies) and these minimums are much more sharply 
defined for curves of small damping ratio. 





FIG. 3 - OUTPUT 0(t) SHOWING VARIABLE INITIAL CONDITIONS 



FIG. 4 - y(t) FOR DELAY IN THE INPUT PULSE SEQUENCE 
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FIG. 9 - {* s f(ZMIN) 
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FIG. 14 - b s f (X) 
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APPENDIX B 


COMPUTER PROGRAMS 


In the programs which follow, the variable X represents 
the Z of the analysis and the variable F represents the I of the 
analysis . 
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